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Abstract 



The leading soft gluon pj- distribution in heavy ion collisions was 
^ I obtained by Kovner, McLerran, and Weigert after solving classical 

^T) ' Yang-Mills equations. I show explicitly this result can be understood 

,—1- . in terms of conventional QCD perturbation theory. I also demon- 

CN i strate that the key logarithm in their result represents the logarithm 

^^ ' in DGLAP evolution equations. 

p^; 1 Introduction 

P-i 

D I In ultra-relativistic heavy ion collisions, after initial collisions, the parton 

^, [ system will be dominated by gluons [jl|. Understanding the distribution of 

soft gluons formed in the initial stage of the collision is particularly interesting 
^ ■ and important for studying the formation of quark-gluon plasma. In terms 

of conventional QCD perturbation theory, a calculable cross section in high 
energy hadronic collisions is factorized into a single collision between two 
partons multiplied by a probability to find these two partons of momentum 
fractions xi and X2, respectively, from two incoming hadrons. The probability 
is then factorized into a product of two parton distributions 0(xi) and 4'{x2), 
which are probabilities to find these two partons from the respective hadrons 
[@]. Because of extremely large number of soft gluons in heavy ion beams, it 
is natural to go beyond the factorized single-scattering formalism to include 
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any possible multiple scattering, and long range correlations between soft 
gluons from two incoming ions. 

McLerran and Venugopalan (MV) developed a new formalism for cal- 
culation of the soft gluon distribution for very large nuclei 0, Q. In this 
approach, the valence quarks in the nuclei are treated as the classical source 
of the color charges. They argued that the valence quark recoil can be ig- 
nored in the limit when the gluons emitted are soft. In addition, because 
of the Lorentz contraction, the color charge of the valence quarks is treated 
approximately as an infinitely thin sheet of color charge along the light cone. 
With these assumptions, the gluon distribution function for very large nuclei 
may be obtained by solving the classical Yang-Mills Equation 0, §. Using 
the classical glue field generated by a single nucleus obtained in the MV 
formalism as the basic input, Kovner, McLerran, and Weigert (KMW) com- 
puted the soft gluon production in a collision of two ultra-relativistic heavy 
nuclei by solving the classical Yang-Mills equations with the iteration to the 
second order [^. The two nuclei are treated as the infinitely thin sheets of 
the classical color charges moving at the speed of light in the positive and the 
negative z directions, respectively. Following this approach, the distribution 
of soft gluons at the rapidity y and the transverse momentum p^ in nuclear 
collisions can be expressed as §, |^ 

dyd-pr (2vr)42ivJ '"^^'"^ ^ p't Wutoff ) - ^^ 

where g is the strong coupling constant, A^c = 3 is the number of the color, 
and Acutoff is a cutoff mass scale 0. In Eq. (|1|), Nq/(2Nc) = SxlJi^, where /i^ 
is the averaged color charge squared per unit area of the valence quark, and 
St is the transverse area of the nuclei. Nq is the total number of the quarks in 
the color charge source. Eq. (||)) is potentially very useful in estimating the 
production of mini-jet rates, and the formation of the possible quark-gluon 
plasma at RHIC [|. 

In the following, I show that this result can be understood in terms of 
QCD perturbation theory. I will also explore under what kind of approxima- 
tion this result matches the conventional perturbative calculation [0]. 




Figure 1: Square of the leading Feynman diagram to the process: qq —^ qqg. 

2 Partonic process qq -^ qqg 

KMW's derivation for Eq. (||) is based on the following physical picture: in 
ultra-relativistic heavy ion collisions, gluons are produced by the fields of two 
strongly Lorentz contracted color charge sources, which are effectively equal 
to the valence quarks of two incoming ions. In order to understand KMW's 
result in the language of perturbative QCD, let us consider a specific partonic 
process: qq -^ qqg, as sketched in Fig. |l|. If we assume that the incoming 
quarks qq are the valence quarks in the initial color charge sources, the par- 
tonic subprocess in Fig. [^ mimics the physical picture adopted in KMW's 
derivation. In this section, I show how to extract the leading contribution to 
the soft gluon production from this diagram. However, as a Feynman dia- 
gram in QCD perturbation theory, the single diagram shown in Fig. |I] is not 
gauge invariant. In next section, I show that with a proper choice of gauge 
and certain approximation, other diagrams are suppressed in soft gluon limit. 

In general, the invariant cross section of the gluon production, shown in 
Fig. |1|, can be written as 

da^g^g/dyd^PT = — \W dps , (2) 

where s = (/i + ^2)^5 and |Mp is matrix element square with the initial-spin 
averaged and the final-spin summed, li and I2 as the momenta of the two 
incoming quarks respectively. In Eq. (^), the phase space 

dps ex -^^{2ir)5{{h-hf)y<j^^{27^)5{{h-k,f) 

x{27r)''5\h + k2-p) , (3) 
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Figure 2: Feynman diagram for 
quark splitting: Pi-^^kiix , kiT < 
Pt)- 



Figure 3: Leading Feynman dia- 
gram contributing to the hard scat- 
tering part H{xli,l2,p)- 



where fci and ^2 are the momenta for the two gluons emitted from the initial 
quarks. Because of the gluon propagators, as shown in Fig. |l|, the matrix 
element square |Mp has the following pole structure: 
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When integrating over the phase space, we see that the leading contribution 
comes from the terms with fc^ — i> or fcg — > limit. Note that k\ and k\ can 
not be zero at the same time, because k\ and ^2 come from different direc- 
tions, and we have the on-shell condition of p^ = (fci -|- ^2)^ = 0. Therefore, 
to calculate the leading contribution, we can first calculate the diagram in 
k\ ^ ^ limit. The total leading contribution is just twice of it, because the 
diagram is symmetric for k\ and /c2. 

To derive the leading contribution at A;^ — ^ limit, we perform the 
coUinear approximation k^ ~ xl\ + OikxT), with k^T ~ ^cutoff « Pt-, 



where Acutoff is a coUinear cutoff scale |]IOl . This approximation means that 
the leading contribution is from the phase space where almost all transverse 
momentum of the final-state gluon comes from the gluon of k2, and fci is 
almost coUinear to h. After such coUinear approximation, the cross section 
in Eq. (^ can be approximately written in a factorized form [§]: 
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where the overall factor of 2 is due to the fact that the leading contribution 
come from two regions corresponding to kf ^ and ^2 ~^ O5 respectively. In 
Eq. (^, Pi^^ki{x,kiT < Pt) represents the probability of finding an almost 
collinear gluon with the momentum fraction x from an incoming quark of 
momentum li, and is represented by Fig. |[ H{xli,l2,p) in Eq- (iD is effec- 
tively the hard scattering between the gluon of fci = xli and the incoming 
quark of I2, and it is represented by the diagram shown in Fig. |[ 

Under the collinear expansion ki ^ xli, the gluon line which connects 
the partonic parts Pi^^ki and H{xli,l2,p) is effectively on the mass-shell, 
and therefore, the partonic parts, P and H in Eq. (|^), are separately gauge 
invariant. The quark splitting function Pi-^^^kiix, fcir < Pt) can be calculated 
in n ■ A = gauge. We have 

Pi^-,kiix,kiT <Pt) 

In n ■ A = gauge, we have the hard scattering function |^ 

1\ /j9+ V 1 1 ^, 2j9+/_ , 



X {xs -2xl^pJf + {2xl^pJ){2p^pJ) . (7) 

For soft gluon production, we define the soft gluon limit as 



f^«l and f±«l. (8) 

Combining Eq. (^, Eq. (|^) and Eq. (^, and taking the soft gluon limit, we 
have [^ 

Eq. @ reproduced Eq. ([l|), which is obtained by solving classical Yang-Mills 
equations, at A'"g = 1. If we consider the total number of the quarks in the 
charge sources of both sides, we need to multiply N"^ to Eq. (Rf). 
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Figure 4: The rest of Feynman diagrams contributing to the process: qq -^ 
qqg, in addition to the diagram in Fig. 1. 

3 Gauge Invariance 

For the gluon production in the process qq — >■ qqg, the single diagram shown 
in Fig. |I| is not gauge invariant. In general, we also need to consider the 
radiation diagrams shown in Fig. ^. Similar to the diagram in Fig. [^, the 
contribution of Fig. ^ and Fig. ^ can also be written in the same factorized 
form: 

(^o-rad ^ I f dx 
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cutoff 



+0(^2|^), (10) 

Pt 

with i = a,b. Pi^^kA^,kiT < pr) is the quark splitting function given by 
Eq. d^). Ha{xli,l2,p) and Hj,{xli,l2,p) are the hard scattering parts from 
the diagrams in Fig. |^a and Fig. |^, and they are represented by Fig. ^ and 
Fig. ^D, respectively. With the contribution from diagrams in Fig. ^, Eq. (|^) 
changes to 

X [H{xlij2,P) + Ha{xli,l2,P) + Hb{xli,l2,P) 
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Figure 5: Feynman diagrams contributing to Ha{xli,l2,p) (a), and 
Hi,{xli,l2,p) (b). 

+interference terms] + 0{ '^^"■'■' ^ , (11) 

Pt 

with the approximation ki = xli + 0{kiT). Feynman diagrams shown in 
Fig. ^ and Fig. ^ form a gauge invariant subset for calculating the hard scat- 
tering parts, H{xli,l2,pys in Eq- (P^- With our choice of gauge, Hi/H ~ 
P-/1- << 1 in the soft gluon limit. Therefore, the contribution from di- 
agrams in Fig. ^ and Fig. ^ can be neglected in comparison with the 
contribution from the diagram in Fig. |I| at fc^ ~ 0. 

Similarly, the contributions from diagrams shown in Fig. ^ and Fig. |d 
can be neglected in the soft gluon limit, when compared with the contribution 
from the diagram in Fig. |l] at k^ ~ 0. Therefore, with the approximation 
of A;^ ~ (or k"^ ~ 0) and the soft gluon limit, and a proper choice of the 
gauge, the contributions extracted from the diagram in Fig. |I| to the leading 
soft gluon production in Eq. (P is gauge invariant. 



4 Discussion 

When we consider the collision between two nuclei, we can treat the two 
incoming quarks in Fig. |I| as coming from two nuclei respectively. In this 
picture, the number of the valence quark Nq is replaced by the quark distri- 
bution in the nuclei. In terms of the parton model, the cross section between 
the two nuclei A and B can be expressed in the following form: 

f?g|ii^/..,.../„.(.0/,/B(..)^^. (12) 

Here zi and Z2 are the momentum fractions of the quarks, and fq/A{zi) and 
fq/B{.Z2) are the quark distributions (or quark number densities) of the two 



nuclei. If we denote pa and pb as the momenta for the two nuclei respectively, 
then zi = h/pA and Z2 = h/pB- Substituting Eq. @ into Eq. (|l5), we have 
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(13) 



where the overall factor 2 in Eq. (^ is now represented by the two terms, 
and S = {pa + Pb)'^ ~ 2pA ■ Pb- In Eq- (|il), ^i = ki/k with i = 1,2, and 
/i = ZiPa and ^2 = Z2Pb- If "we denote the momentum fraction of gluon ki 
with respect to pa as z[ = ki/pA, and ^2 with respect to ps as Zg = ^2/^5, 
we can rewrite Eq. (^) in terms of z[ and ^2- 
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(14) 



According to the QCD factorization theorem 0, we see that the part inside 
the square brackets is actually the gluon distribution from nuclei A (or B) at 
the factorization scale /i|. = p"^, with only the quark splitting function Jill, 

fg/Aiz[,iJ,l = Pt) = —^ fg/Aizi)Pi,-,kiiz[/zi, kiT < Pt) 

J Zi 



+term from gluon splitting. 



(15) 



Using Eq. (JTSj), we can then reexpress Eq. (JT^) as: 

dydpl ^ 2(27r)3 25 
dz' dz 



X 



— [fg/Aiz',fJ.l = PT)fg/B{z)H{z'pA,ZpB,p) 

J A/ A/ 

+fq/A{z)fg/B{z',fJ^l = pI)H{zpa,z'pb,p)] , (16) 

which is the factorized formula for two-to-two subprocesses in the conven- 
tional perturbative QCD for the nucleus-nucleus collisions. In KMW for- 
malism, only the valence quark color charge was used as the source of the 



classical charge of colors. As a result, the gluon splitting term in Eq. (15) is 
neglected for the distribution fg/A- 

From the above comparison, we conclude that by solving the classical 
Yang-Mills Equation to the second order in iteration, KMW's result is con- 
sistent with conventional perturbative QCD at the leading logarithmic ap- 
proximation. The logarithmic dependence shown in KMW's result basically 
describes the logarithmic DGLAP evolution of the quark distributions [|T^]. 
However, the iteration in this approach is different from the expansion series 
in conventional perturbative QCD. The McLerran-Venugopalan formalism 
was later further developed to include the harder gluons into the charge den- 
sity /i^ and treat the charge source as an extended distribution which depends 



on the rapidity |jT3[ . It will be potentially very useful if this new approach, 
after including higher orders of iteration, can include the parton recombina- 
tion [|T^ and other non-perturbative effects which are not apparent in the 
normal perturbative calculation. 
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